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Abstract 

In this article, we study the problem of designing a Decoherence Control for quantum systems with 
the help of a scalable ancillary quantum control and techniques from geometric control theory, in order 
to successfully and completely decouple an open quantum system from its environment. We re-formulate 
the problem of decoherence control as a disturbance rejection scheme which also leads us to the idea of 
Internal Model Principle for quantum control systems which is first of its kind in the literature. 

It is shown that decoupling a quantum disturbance from an open quantum system, is possible only with 
the help of a quantum controller which takes into account the model of the environmental interaction. 
This is demonstrated for a simple 2-qubit system wherein the effects of decoherence are completely 
eliminated. The theory provides conditions to be imposed on the controller to ensure perfect decoupling. 
Hence the problem of decoherence control naturally gives rise to the quantum internal model principle 
which relates the disturbance rejecting control to the model of the environmental interaction. 

Classical internal model principle and disturbance decoupling focus on different aspects viz. perfect 
output tracking and complete decoupling of output from external disturbances respectively. However for 
quantum systems, the two problems come together and merge in order to produce an effective platform 
for decoherence control. In this article we introduce a seminal connection between disturbance decoupling 
and the corresponding analog for internal model principle for quantum systems. 

I. Introduction 

Quantum Information and Quantum Computation hold the key to faster information processing and 
better and reliable communication (T). The properties, the quantum superposition, coherence and en- 
tanglement are vital to quantum information processing. Quantum measurements in general collapses a 
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quantum state into set of bases decided by the observable. Decoherence (2l is the process by which the 
quantum system loses the coherence and superposition by continually interacting with the environment. 
A quantum system that is continuously interacting with the environment is called an Open Quantum 
System 0. Decoherence is conceptually equivalent to a continuous and forcible collapse of the wave 
function of the system onto the basis decided by the environment(also called the pointer basis (H). In 
practice, this adiabatic process takes a finite time in the order of a few milliseconds thus rendering 
the quantum system classical. The problem of decoherence is currently the biggest roadblock towards 
exploitation of quantum speedup in computation. Thus far, many researchers have proposed multitude of 
ways to control decoherence in such open quantum systems, of which a few representative contributions 
include open-loop pulses @ @ O, and control within Decoherence Free Subspaces J9). Open loop 
pulses techniques are designed for systems that are acted upon by pre-programmed control pulses. Such 
methods also have the tendency to annihilate useful non-predetermined control in addition to suppressing 
decoherence. This is debilitating for quantum information processing and do not work under arbitrary 
and undetermined control. Another class of ideas is based on the Decoherence Free Subspace(DFS) 
which are proven immune to decoherence due to the degeneracy of the basis vectors with respect to the 
decoherence interaction. Such methods aim at encoding and steering the quantum information within such 
a subspace at all time. Such a strategy does not admit arbitrary control Hamiltonians as any transition 
out of the subspace would subject the state to decoherence and hence loss of information. 

Another class of work which is based on symptom or syndrome correction is error correction codes 
|[T0l ifTTI . These methods aim at correcting the observed effects or symptoms of decoherence. These 
methods usually require ancillary bits to encode a specific quantum information in a redundant fashion 
and perform posterior transformations depending on the observed error syndrome. Such methods require 
number of ancillary/redundant bits proportional to the size of the original system and might not be 
scalable in the long run. 

Hence most of the proposed techniques are either ad-hoc, or limited in control functions or not scalable. 
In this work, we propose a scalable, strategy which preserves action of useful controls and lets the system 
evolve according to the same while eliminating the effects of decoherence. This is applicable to a wide 
class of control as well as decoherence Hamiltonians. 

Our work is orthogonal to almost all published work in the literature on decoherence control/quantum 
disturbance decoupling. We strongly believe that our study provides another avenue of research for 
decoherence control, informs the readers in this field of new directions to be investigated for the same 
problem and lends a deeper insight to quantum disturbance decoupling. Needless to say, reviewing the 
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extensive nature of the work in the literature on using density matrix approach to study the problem, 
one cannot expect to give a complete solution to this important problem in one paper. We present some 
pivotal and important results on decoherence control in this paper with the following main contributions: 

1) Open Loop Invariance: Utilize differential geometric tools to perform structural analysis and 
extract important information regarding susceptibility of the given system to decoherence. This 
helps determine a priori whether or not the given system is immune to decoherence interaction 
and can be used to avoid going through, tedious, sometimes futile and time consuming work. 

2) Active Controller: Provide results in terms of the given control equation and available control 
resources, whether complete decoupling of the effects of disturbances is possible, with the help of 
an active controller. 

3) Ancillary Quantum Control: If the system is not decouplable, design an effective control system 
via an ancillary quantum system and an active controller, that achieves complete decoupling under 
arbitrary and non-predetermined control. To the best of our knowledge the decoherence prevention 
in the presence of arbitrary useful control has not been addressed before. Moreover our results also 
provide a systematic way(not ad-hoc) to construct the desired control. 

To this effect, we first provide the criteria for any system to be naturally immune to decoherence in terms of 
the Lie Algebra of the operators involved, in the presence of arbitrary user generated control. The treatment 
is powerful and general enough to yield Decoherence Free Subspace(DFS) as a special case of the open 
loop control. In addition, this yields best ways to encode a given quantum information that is immune 
to decoherence under arbitrary control. Secondly, for those systems that are not immune to decoherence 
under arbitrary control, and systems undergoing decoherence, we employ an active controller. At this 
point, we transition from an operator algebra method to a vector field method on the tangent space of the 
manifold as this offers additional valuable insights into the geometric nature of the problem. We present 
a scalable construction involving an ancillary system(single ancillary qubit for a finite number of system 
qubits) to achieve complete decoherence control. All of the analyses mentioned above are performed in 
the presence of arbitrary user generated control which preserves useful work while eliminating only the 
effects of decoherence. 

Finally, we present the simulation results with the above control strategy. The above mentioned ideas 
come together in a coherent way into what is called "Quantum Internal Model principle" wherein the 
model of interaction with the environment is indispensable to efficient disturbance decoupling which will 
be discussed in the last section of the paper. 
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II. Mathematical Preliminaries 

A pioneering effort to study quantum control systems using bilinear input affme model was carried 
out by Huang et. al |[T2l . The model has since found various applications and is extremely useful in 
analyzing the controllability properties of a quantum system on the state space of analytic manifolds 
|[T3l which draws upon the previous results on controllability of finite dimensional classical systems by 
Sussman and Jurdjevic |[T4l which in turn follows the results by Kunita lfl31 |[T6l and Chow ifTTl . In this 
section we explore the conditions for a scalar function represented by a quadratic form to be invariant 
under the dynamics of the above model(with the additional assumption of time-varying vector fields) in 
the presence of a perturbation or interaction Hamiltonian. Such a formalism can be seen to readily relate 
to decoherence in open quantum systems wherein a perturbation Hamiltonian that couples the system to 
the environment can be seen to play the role of disturbance. Classical disturbance decoupling |[T8l |[T9l 
|20l provides insightful results on eliminating the effects of disturbance from output, however it will also 
be seen that the aforementioned is not quite similar to quantum decoherence control problem and one 
should be extremely careful in adapting the classical results to quantum regime. 

Let the quantum control system corresponding to an open quantum system [3l interacting with the 
environment (figure (Q])) be given by, 
df(t x) 

K m ' = [H ® Ze(t, x) + l s <8> H e (t, x) + H SE (t, x) + £l=i Ui{t)Hi ® l e (t, x)]£(t, x) (1) 

where, T-L s is the system's Hilbert space and 7i e the environment's Hilbert space. T~L S could be finite or 
infinite dimensional and % e is generally infinite dimensional. x) is the wave function of the system 
and environment. Hq and H e are operators corresponding to the drift Hamiltonian of the system and 
environment while Hi's correspond to the control Hamiltonian of the system. H$e governs the interaction 
between the system and the environment. The above operators are skew Hermitian and assumed to be 
time varying and dependent on the spatial variable. In addition to the above dynamical equation, we 
introduce a complex scalar functional, y(t), as a bilinear form that carries information about the system, 

y(t,Q = (£(t,x)\C(t,x)\S(t,x)) (2) 

The bilinear form of the function resembles the expected value of an operator, but it does not necessarily 
correspond to a measurement output. The operator C(t, x) (Hermitian or non-Hermitian) is a time-varying 
operator acting on system Hilbert space. In all the subsequent analysis we study the invariance properties 
of the above scalar map of the system acted upon by decoherence interaction. This could be thought of 
as a function that has to be regulated in the presence of controls and disturbance. 
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Fig. 1. An Open quantum system interacting with the environment via Hse 

Definition II.l. Let y(t,£) = f(t,x,ui,--- ,u r ,HsE) be a complex scalar map of the system as a 
function of the control functions and interaction Hamiltonian over a time interval to < t < t\. The 
function is said to be invariant of the interaction Hamiltonian if 



for all admissible control functions u\, ■ ■■ ,u r and a given interaction Hamiltonian Hse- 

Forms of y(t). The above equation takes a quadratic form in the state |£) of the combined system 
and the environment. Some of the possible functions of interest, chosen for regulation are, 

(i) An expected value of a physical observable or an observation. The operator C, if Hermitian, can 
be a non-demolition observable and y(t) is the output of the measurement performed on the system. 

(ii) By a suitable choice of the operator C the value y(t) can represent the coherence between the 
basis states of interest. It can be seen that a suitable value of the operator C could yield the off-diagonal 
terms of the density matrix for the function y(t). 

For example C = \si)(sj\ g) I e can be seen to yield the coherence between the orthogonal states 
of the system |sj) and \sj). For the pure state £ = ^ C «N)> v(t) = c l c j anc ^ f° r tne entangled state 
£ = ^2ci\si)\ei) where | e-i) are the orthogonal states of the environment, a similar calculation yields 
y = 0, denoting no coherence between the basis states. 

{in) The operator C could also be a linear combination of multiple non-Hermitian operators denoting 
coherence information of multiple basis vectors. This form for the operator C is extremely useful in 
studying invariance of quantum information under external influences. This is presented in detail in the 
section ^ on DFS. 



f(t,X,U!,--- 



u t ,Hse) = f(t,x,ui, ■ ■ ■ 



u r , 0) 



(3) 
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Let us define the corresponding free, control and the interaction vector fields as follows. 

K =(H + H e )\C) (4) 

Ki =Hi\C) (5) 

Ki =H SE \0 (6) 

Here we have suppressed the dependence of the Hamiltonians on time and spatial variable. As most of the 
practical systems are time-invariant and locality-invariant, this is a reasonable assumption. The following 
lemma ET1 provides the basic conditions necessary for invariance of the scalar map with respect to the 
interaction vector field. 

Lemma ILL Given that the quantum control system ([7]) is analytic on the analytic manifold, the 
corresponding scalar map given by equation (|2]) is invariant under given HsE(or the corresponding 
vector field, Kj ) if and only if for all integers p > and any choice of vector fields X± , • • • , X p in the 
set {Kq,Ki,--- ,K r }, 

L Kl L Xl ■ ■ ■ L Xp y{t, = 0;for all t, £ (7) 

Lemma III. 1 1 implies that the necessary and sufficient conditions for the scalar map y of an analytic 
system to be invariant of the interaction vector field, Kj are, 

L Kl y(t,£)=0 

L Kl L Kio ---L Kin y(t,O = (8) 

for all t, £, for < 

*0: ■ ■ ' ) in 5: r an d fi > 0, where Kq, • • • , K r are the corresponding drift and control 
vector fields and Kj, the interaction vector field. 

III. Invariance for the Quantum System 

With the preceding mathematical preliminaries in place we can now apply the above conditions to 
the quantum system with careful consideration to the nature of the complex functional and the analytic 
manifold. In this section we present two different cases, (a) the open loop invariance, (b) open loop 
invariance under an external controller. Both these cases are studied for an open quantum system acted 
upon by arbitrary useful control. 
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A. Open Loop Invariance with Arbitrary Control. 

We can now state the condition for invariance of the scalar function y(t) with respect to a perturbation 
or interaction Hamiltonian, the proof and motivation for which is presented in ETTl . 

Theorem III.l. Let Co = C(t) and for n = 1, 2, • ■ • , define 

C n = span{ad J H C n -i(t)\j = 0, 1, . . . ; i = 1, . . . , r} 

Cn = l (adH + ^j C n ;j = 0,l,-- - 1 



where H = H$ + H e , the drift Hamiltonian of the combined system and environment, and Hi, i = 1 • • • , r, 
the control Hamiltonians. Define a distribution of operators, C(t) = span{Ci(t),C2(t), ■ ■ ■ ,C n (t), ■ ■ ■ }. 
The scalar function represented by equation (f2]) of the quantum system is decoupled from the environmental 
interactions if and only if, 

[C(t),H SE (t)} = (9) 

Proof: The proof follows by noting the equivalence of equation ([8]) with the above condition. Consider 
the following term Lx zg ■ ■ ■ Lx ik y(x) for any k > 1, and io, ■ ■ ■ G {0, ■ ■ ■ r}. It can be expressed as a 
bilienar form of an operator of Lie brackets of Hi ,Hi x , • • • Hi r , C and their time derivatives as follows. 
In particular for k = 0, 

L Kio y = {Z\[C,H io ] + 5(i )jC\0 4 (flr^) 

where 5(io) is the delta function that takes value 1 when zq = and the operator Ti as defined above, 
is such that T\ G C±. Similarly for k = 1 we have 

Lk^Lr^U = (£\llC,H l0 ],H tl ] + Mio^CtfiJ +*(ti)^([C J iZi ] + «5(io)-|c)|0 

= mio 

and T2 € C2. Continuing so, in general we have T n G C n . Now via condition (HJ, we have (£| T" n ] |£) = 
0,V£, or [i^5B)^n] = in general for invariance. Since the condition is true for any n > and any 
T n and since the vector space of bounded linear operators is complete we have [Hsei a i^i\ = 
a i[HsE, Ti] = for Oj € JR. The converse is true by noting that any operator in the distribution 
C (i.e) for any T G C can be decomposed into a sum of operators ^ a{Ti for Tj G Ci and given 
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[Hse, YHHLo a iTi] = OVoj which is true only when [HsE,T n ] = for any n. Hence from the previous 
equations L Kl L Kln L Kln _ 1 ■ ■ ■ L Kzo = for i , ■ ■ ■ , i k € {0, • • • r}. ■ 
We now present two qualitatively different examples, a system undergoing decoherence and a system 
that is immune to decoherence due to its Decoherence Free Subspace, to illustrate the applicability of 
the above open loop invariance theorem to practical quantum control systems. 

1 ) Electro-optic Amplitude Modulation: Consider a driven electromagnetic system in a single mode 
subject to decoherence. The control system describing the oscillator under semiclassical approximation 
is 

j t £(t) =(utJa + E Ujbjbj + iu(t){a) -a) + aJ2 9jbj WE tofyft) 
j j j 

where the system represented by mode a is coupled to a bath of infinite number of oscillators, bj 

with corresponding coupling constants gj and where £(t) is the combined wave function of the system 

and bath. Here a, bj and cb^bj, denote the photon creation and annihilation operators respectively. The 

control u(t) is the strength of the input current and Hq = ooa^a + J2j u jb]bj, Hi = {a) — a) and 

Hse = a Ylij 9jbj + ^ Ylj 9jbj are the drift, control and decoherence Hamiltonians respectively. Let the 

system be monitored by a non-demolition observable E2l . 

C(t) = aexp(iwi) + exp(— iut) 

with the corresponding output of the non-demolition measurement given by y(t) = (£(t)|C(t)|£(i)). 
Following theorem IlII. 1 1 we have [C(t),H\] = e lult + e~ tLOt = 2cos(wi) with vanishing higher order 
commutators. Hence C\ = {ciC + C2lcos(u>t), Vci, C2 G R} and since [C(t), Ho] + dC/dt = we have 
Ci = C\. Since the commutant of the interaction Hamiltonian Hse = a Ylj 9jbj + a ^ Ylj 9j"j with the 
elements of the set Ci not all zero, the non-demolition measurement is (i) not invariant of the interaction 
Hamiltonian, (ii) no longer back action evading due to the presence of the interaction. The measurement 
of the observable C(t) would thus reveal information about the decoherence of the system. 

2) Decoherence Free Subspaces(DFS) of a collection of 2-level systems:: The above theorem 
can also be applied to the problem of analyzing the DFS 0. Decoherence free subspaces camouflage 
themselves so as to be undetected by the interaction Hamiltonian due to degeneracy of their basis states 
with respect to Hse and the special algebraic properties of the interaction Hamiltonians. For a collection 
of 2-level systems interacting with a bath of oscillators the corresponding Hamiltonian is, 

N N 

h = f E ^ + E "4** + E E WI + 9th) 

j=l k k j=l 
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where the system is assumed to interact through the collective operator Y^j a z and g^s describe coupling 
to the mode k. An inquiry into what information about the system is preserved in the presence of the 
interaction can be answered by expressing the operator C acting on the system Hilbert space in its general 
form in terms of the basis projection operators, 

c(t)= ^01 

i,j=0..2 N -l 

and solving for condition ©. For a simple N=2 system, we have after straight forward calculations 
C = span{^ c l3 \i)(j\.(j^ - i« + j(2) _ iWfy K = 0, 1, 2...} 

where j'W denotes the I th letter (0 or 1) of the binary word j. Equation © which is, [C,Hse] = 
implies, 

Y^<H j \i){j\.U il) -i {1) +j (2) -i i2) ) K = 0,yK = 1,2,3... 

or non-trivially j'W = +i^ 2 \ or the two words have equal number of l's. The above calculations 
are valid for any finite N, a specific example for N = 3 is C = |000)(000| + |001)(001j + |010)(100| + 
|011) (101 1. Of particular interest are terms such as, |011)(101| and |010}(100| as the corresponding 
y(t) = (C(t)\C(t)\^(t)} which is a function of the coherence between the basis states |011) , 1 101) and 
1 010) , 1 100) is predicted to be invariant under the interaction. It is worth noting that the operator C(t) 
acting on system Hilbert space here need not necessarily be Hermitian and only describes the quantum 
information that is preserved. 

(i) 

Decoherence in the presence of control: In the presence of the external controls Hi = UiU x , the 
invariance condition (© is no longer satisfied for the operator C = ^ J=0 2 N -i K)0M 7^ 3 as 
[[C, a^],ai j) ] 7^ and hence the coherence between the states \i), \ j) is not preserved. This is because 
of the transitions outside the DFS caused by the control Hamiltonian. The above formalism is helpful 
in analyzing in general, class of information that would be preserved in the presence of interaction 
Hamiltonian which in turn could tell us about how to store information reliably in a quantum register 
in the presence of decoherence. Hence, in contrast to passive decoherence avoidance in the absence 
of external controls, this approach can be used to determine the prudent means to encode quantum 
information, that stays immune to the decoherence interaction, even in the presence of arbitrary controls. 

In summary, the notion of open loop invariance, (a) naturally gives rise to DFS described by the 
operator C, and helps perform extended analysis on the same, (b) was used to determine if a given scalar 
function(in this case non-demolition measurement) was affected by decoherence, and (c) could be used 
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to design an operator C, which under the given system Hamiltonians and decoherence interaction would 
generate an invariant scalar map. 

B. Controller: u = a(£) + /3(£,)v 

In this section we consider the case where the scalar map is not invariant in the presence of decoherence 
control(for eg. the coherence between states 1 10) and |01) of a 2 qubit system) and investigate the role 
of an active controller in order to decouple the same. From this point onward(and all the subsequent 
sections) we assume that the scalar map is not invariant, and that the decoherence operator(i7s e) affects 
the scalar map. We study the effectiveness of an external controller whose control is dependent on the 
state of the system, in order to achieve invariance of the scalar map. The controller is assumed to be 
of the form, u = a(£) + (3(£)v where the u is implemented as a transformation involving matrices a 
and p. Here v is 1 x r, u is 1 x r, a(£), is a 1 x r vector and /?(£) is a non-singular matrix of of size 
r x r(where r is the number of open-loop controls). 

The realizability of the controller and implementation of the obtained control is a part of the ongoing 
work. It could potentially be realized with the help of a quantum machine(coherent control) or via 
quantum measurement/estimation theory. As outlined earlier, the focus of the current work is geometric 
analysis and a method for decoherence control design. The realization of the controller is an important 
open problem currently under study which could yield a few dissertations by itself. 

The above form of an external control(u = a(£) + f3(£)v) is general enough to encompass all popular 
control strategies as well as preserve the input-affine form of the original quantum control system even 
after the application of control. Consider the following system that is acted upon by the above controller 
of the form u = a(£) + /3(£)v, 

d r r 

x) = {H + H e + J2 <XiHi)£(t) +J2 Vi J2 &i H £<t) + HseZ® (10) 

i=l j=l 

where the new drift vector field, K = (H +H e +^2 ctiHi)£(t), control vector fields = ^ ■ /3ijHj£(t), 
and decoherence interaction Kj = HsE£{t), are identified for the controlled system. 

As stated earlier, the necessary and sufficient conditions for a scalar function y(t) of the system to be 
invariant of the interaction vector field are, 

L Kl y(t) = 

L Kl L kio ---L ki j{t) = Q (11) 
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for < 

^0) ■ ' ' i *n — r an( i n — 0- The above conditions can be stated in the form of the following result, 
which states it in terms of the operators defined in theorem dIII.ll ). 

Definition III.l. The scalar map (0 of the system (0) is said to be decouplable if there exist control 
parameters a and (3, such that, under the corresponding controlled system < \10i , the scalar map ((2]) is 
invariant in sense of definition Ml. 1 D . 

We now state the condition for the decouplability. 

Lemma III.2. For the scalar map ([2]) of the quantum system ([7]), is decouplable only if, 

[C,H SE }=0 

[C(t),H SE (t)]cC(t) 
where the distribution C(t) is as defined in theorem( J I ). 

Proof: The above lemma only provides the necessary condition for the invariance. In the equations 
below we suppress the summation symbol and follow Einstein's convention, wherein a summation has 
to be assumed whenever a pair of the same index appears. Expanding out the corresponding terms of the 
equations (fTTT ). where the following equalities must hold for all £, 

L Kl y = {S\[C,H SB ]\0=0 

L Kl L ki y = (£\[[C,f3ijHj],HsE] + [C,Hj]L K Aj\0 = 

L R L Ro y = UCJuHi] + [[C,H + n 1,,!!, + [C, Hj]L R atf) = 

L i<i L R z L k y 

= mC, PuHiIHse) + [[C, HAL^, H SE ] + [[[C, H + a.Hj], (3^}, H SE ] 

+ [C, Rr\L Kl P iX + [C, HjlLx^aj + [[C, H^H^Lk^u + [[C, H^H^L^aMO 
=0 (12) 

The last equation above, provides a set of simultaneous equations to solve for the control parameters 
a and /3 of the active controller, in order to achieve invariance. The above equation contains two types 
of terms. The terms containing H$ E and terms that do not. The terms whose commutator with H$ E is 
computed, is found to belong to the distribution [C(t), H$ E ] whereas the terms without Hs E is seen to 
belong to C(t). The above calculation can be extended to finite number of terms to arrive at the result. In 
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order for the above equality to hold, in general one finds that, the condition for decouplability is relaxed 
to, 

[C(t), H SE ] C C(t) (compare to eq. gj (13) 

■ 

However, in order to solve equation (fT2l ) and consequently equation (TT3T > for the control parameters(a 
and f$), it is important to study the properties of the operators in C. It can be seen that the distribution C, 
is generated by operators acting only on system Hilbert Space (C, Hq, Hi - - - H r ), whereas the operator 
H$e acts on the joint, system + environment Hilbert Space. Therefore, the above equation ( fT3l cannot 
be solved for the control parameters unless [C(t), H$e] = 0. This is same as the open loop invariance 
without the active controller, which implies that as long as H$e and Hq, ■ ■ ■ , H r , act on different Hilbert 
spaces, the controller cannot act as an effective tool in decoupling the system. 

Alternatively, for the controller to be an effective tool in solving the decoherence problem, the control 
Hamiltonians Hi 's have to act non-trivially on both the Hilbert spaces which would enable all the operators 
in equation ( fT3l to act on the joint system-environment Hilbert space. 

In the rest of the paper we will outline a construction, involving an ancillary system and the active 
controller of the form u = a(£) + /3(£)i> in order to decouple and achieve complete invariance. We will 
revisit the 1 and 2-qubit systems and present the applicability of the construction in achieving the final 
goal of decoherence control. 

The operator algebra method outlined above was helpful in arriving at the invariance condition for 
open loop, with and without the active controller. However, it only provides the necessary condition in 
order to be able to achieve invariance under the action of the controller. Hence at this point we resort to 
an alternative approach to analyze the same problem, via the invariant subspace within the tangent space 
of the analytic manifold. In contrast to the operator algebra approach, which was based on the operators 
Hq, Hi,-- - ,H r and H$e, we now use the control vector fields K\, ■ ■ ■ , K r , and the decoherence vector 
field Kj, in order to analyze the invariance of the function, y and leverage the geometry of system on 
the analytic manifold. 

IV. Invariant Subspace Formalism 
In this section we present a alternate formalism to analyze the invariance of the function y. 

Definition IV. 1. Any vector field K T = Kj satisfying equations ([§]) is said to be in the orthogonal 
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subspace of the observation space spanned by the one-forms, 

O = span{dy(t, £),dL Kio y(t, £), • • • , dL Kio ■ ■ ■ L Kin y(t, £)>•••} 

VO < «0j ■ " ' j *n < r a?I <^ n > (14) 

Denoted by K T G 0^ 

Lemma IV. 1. 77?e distribution O 1 - is invariant with respect to the vector fields Kq, ■ ■ ■ ,K r under the 
Lie bracket operation, i.e., if K T G 0^~, then [K T , Kj\ G O 1 - for i = 0, • • • , r 

Equations © after subtraction imply L^ Lx 7 y(t) — LK I LK y{t) = Ly Ka Kl ^y(t) = 0. Similarly it is 
possible to derive other necessary conditions viz. L\ KotKl ^LK j y{t) = and LK j Ly Ko Kl T i y(t) = for 
invariance ET1 which in turn imply that L^ Ko ^ Ki ^ K .^y(t) = 0. In fact the above pattern of equations can 
be extended to any number of finite Lie brackets to conclude that, 

L[[...[K I ,K H ],K i2 ]...K ik }y(t) = (15) 

1 < ii,*2j "•' ,ik — r > which leads us to the definition of an invariant distribution A of vector fields 
with the following properties, 

K v G A =>■ L^y(i) = (16) 

K„ K^A => K,] G A (17) 

and for any control/drift vector field Kq : K\, ■ ■ ■ , K r , 

K v e A [K v ,Ki] G A,Vt£0,--- ,r (18) 

This distribution is involutive and it is also observed(from the definition) that K[ G A. Such a distribution 
A is contained within ker(dy(t, £)). Hence i£/ G A C ker(dy). From the necessary conditions listed 
above the distribution is invariant under the control and drift vector fields Kq, ■ ■ ■ ,K r . Simply stated, 

[A,Ki] c A,Vt €(),■■• ,r (19) 

It is also to be noted that the above calculations are reversible and the original necessary and sufficient 
conditions ([8]) can be derived starting from the invariant distribution. Hence the necessary and sufficient 
conditions for open loop decouplability can now be restated (without proof) in terms of the invariant 
distribution |[23l . 

Theorem IV.2. The scalar map y(t) is unaffected by the interaction vector field Kj if and only if there 
exists a distribution A with the following properties, 
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Fig. 2. The surface represents equal value of y(t) for different values of |£) and the corresponding vector fields at the point |£) 
is denoted by the arrows and nullspace ker(dy) is the tangent to the surface at the point The dotted arrows lie in ker(dy) 
and the solid arrows need not necessarily be contained in the same. The necessary condition for open loop invariance requires 
that H SE \0 G ker(dy). 

(i) A is invariant under the vector fields Kq,K±, ■ ■ ■ , K r i.e, 



(ii) Ki G A C ker(dy(t)) 

A geometric representation of ker(dy) is illustrated in figure ©. The existence of the invariant subspace 
A is essential to decouplability of the given system. We now analyze the decouplability with the help of 
the active controller. 

A. Active Controller: u = ck(£) + P{£)v 

In this section, we study the synthesis of the control parameters in the fundamental limit, that ensures 
complete decoupling from Hse- 

Definition IV.2. A distribution A is said to be controlled invariant on the analytic manifold if there 
exists a control pair (a, (3), a, vector valued and f3, matrix valued functions such that 



[A,Ki] C A,Vt GO,-- - ,r 



(20) 



[ir ,A](£)cA(£) 



(21) 



[K h A](0 c A(£) 



(22) 



r 



r 



where, Kq = Kq + 
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The vector fields Kq and Ki are the new drift and control vector fields of system under the action of 
the controller (a, /?). The above definition of controlled invariance is a simple extension of the invariance 
condition d20l ) for the open loop case. It is now possible to express the necessary and sufficient conditions 
for the controlled system (Kq, K\, ■ ■ ■ , K r ) to be decoupled from the interaction vector field Kj just as 
we were able to provide conditions for open loop decouplability. The following theorem provides the 
necessary and sufficient conditions. 

Theorem IV.3. The scalar map y(t,£) = (£|C(i)|£) can be decoupled from interaction vector field Kj 
via suitable analytic control parameters (a(£), /?(£)) if and only if there exists an involutive distribution 
A defined on the analytic manifold, such that, 

[4A]CA + G (23) 

[*Q,A]cA + G (24) 

and Kj € A C ker{dy) and where G = span {K\, • • • , K r } 



The proof of the above theorem invokes a construction for the sufficiency 11231 . which also provides the 
means to synthesize the control parameters and /3(£). It is based on this construct that we determine 
the invariant subspace for the given system and utilize the same to synthesize the parameters. We shall 
now examine the decouplability of 1-qubit and 2-qubit systems via the conditions stated above and also 
present an application of the theorem in designing the control for the 2-qubit system in section dXl) . 

V. Examples 

We are interested in studying the decoherence of 1-qubit and 2-qubit systems in the presence of 
arbitrary controls. Although the theory developed here is general enough to be applicable for all scalar 
maps represented as bilinear forms, we focus our attention to the decouplability properties of the scalar 
map which represents the coherence between a set of representative basis states. For the 1-qubit system 
it is given by C = |1)(0| and for the 2-qubit system it is C = |10)(01|. The 1-qubit and 2-qubit system 
coupled to the environment are modeled as Spin-Boson systems. Consider the 1-qubit system, 

^fp = y^e(i) + X) W * 6 I 6 *£(*) + + «2°»£(t) + Yl a ^9kb\ + 9* k b k )C(t) 

k k 

with the scalar map, y(t) = (£(t)|C|£(i)), where C = |1)(0|, the coherence between the states |0) and 
|1). The open loop controls, u\,U2 acting via the Hamiltonians a x , a y are assumed piecewise constant. We 
now examine the necessary and sufficient conditions for decouplability via the two approaches developed 
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so far: Opearator Alegbra: It can be seen that for the 1-qubit open quantum system, [C, Hse] = is not 
satisfied for the given, C = |1)(0| and Hse- This is vital for both open loop invariance under arbitrary 
control([C, Hse] = 0, eq. (0) as well as decouplability via an active controller (lemma HH.2| ).Invariant 
Subspace: The necessary condition for open loop invariance(K/ E ker(dy), threorem (1IV.2I) . cond. (ii)) 
is again not satisfied by the 1-qubit system, as Kj = Yl< k a z{9kb\ + 9 k b k )i{t) ^ ker(dy(t)) because 
LidUit) 7^ 0. This is also vital for decouplability via an active controller (theorem (IIV.3I )). Hence 
the conclusion that an open 1-qubit system is not decouplable is independently arrived at by both the 
formalisms. 

Now, consider the following open 2-qubit system, 

lew) +E [ E°*° I (9kbl+ 9 * k b k )\m) 



dt 1^2 

, j=l k I k 

+ ( Ul (t)aW + u 2 (t)a^ + U3 (t)4 2) + u A (t)of))Ut)) 

with the scalar map, y(t) = (£(t)|C|£(i)) where C = |01)(10| and the given Hse- Operator Algebra: 
As was previously shown in section (IIII-A2I ). this system has a DFS of dimension 2, span{|01), |10)}, the 
states within which remain coherent in the absence of controls. It is also seen that [C, Hse] = 0. While 
this part of the necessary condition is satisfied, the sufficient condition for open loop invariance under 
arbitrary control([C, Hse] = 0, eq. ©) as well as decouplability via an active controller (lemma Hll.21 ) 
is problematic as outlined in section (IIII-A2I ) and the discussion following lemma (IIII.2I ). Invariant 
Subspace: Again, it can be seen that the interaction vector field Ki = Yljk a z G?fcfyfc + 9k^k)^(t) 
belongs to ker(dy(t)), (because LK T y(t) = 0) as required by theorem ( IIV.2I ) for open loop invariance. 
However for necessary and sufficient conditions (l23T)-(f24l). we compute the Lie bracket of the control 
vector fields, K\,K2 £ G, with Kj € kei(dy) C A , and note that, 

[k 112 ,kj] = k«io,E^ ) o?4+ 9p k )\0] 



c-^^l(9kb{ + glb k )\0, (25) 



y 

k 

up to a constant c (where | in the subscript is a placeholder). It can be seen that, [-fT^, -K'/] does not 
belong to ker(dy)(and hence A), nor does it belong to the control distribution G. Hence [i^x^,-^/] does 
not belong to A as required by open loop invariance, nor A + G required for decouplability via an active 
controller. Again the observations made with the help of the two formalisms coincide implying that no 
collection of two level systems are open loop invariant under arbitrary control nor can be decoupled with 
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an active controller in its present form. However, it can be noted from the aforementioned discussions, 
that by suitably changing the control Hamiltonians Hi(or the control vector fields Ki), it is possible to 
satisfy the necessary and sufficient conditions. Although, the coherence operator C and the interaction 
Hamiltonian H$ e are fixed for a given system, we can modify the control vector fields in order to meet 
the necessary and sufficient conditions. It is for this reason that the 1-qubit system is not decouplable(as 
[C, Hse] 7^ 0) but the 2-qubit system can be decoupled via suitable modifications. The rest of the paper is 
devoted to studying the means by which such a modification can lead to complete decoherence control. 
For this purpose, we employ a scalable construction and the quantum controller, the single ancillary 
qubit in order to affect such a transformation. Henceforth, we will confine ourselves to the study of 
decouplability of the 2-qubit system and present the ensuing quantum internal model principle. 

VI. An Ancillary Quantum Controller 

Consider the following construction employing an ancillary qubit as a quantum controller, with the 
additional property that its decoherence or strength of the environmental interaction can be modulated 
externally at will. With this construction it is necessary to maintain only one qubit within a modulated 
environment, as opposed to a whole system of finite number of qubits. This greatly simplifies the 
realizability and operability of a practical quantum computer in an ambient setting, without the need 
for large supercooled environments. This system is now allowed to interact with our qubits of interest 
through an Ising type coupling Ji, J2 (figure ([3])). Such a hardware is currently under investigation 
and development 11241 with encouraging experimental results. The state vector is now the total wave 
function of system+ancillary+environment. Both the qubit systems are assumed to interact with the 
common environment with the only additional requirement that the ancillary qubit's decoherence rate 
be controllable. Physically this amounts to a coherent qubit with controllable environmental interaction. 
The scalability and advantages of this construction are analyzed in the next section. The control system 
governing the mechanics following the Schrodinger equation (1261 ) is given by, 

\j=i k j j,k 

+ u 3 (t)4 2) + «4(t)4 2) + y4 6) + «s4 6) + u ^i b) + urJiaPa^ + u 8 J 2 ai 2 Vf) f (i) 



u 9 



E4 6 Wk + «£W*) ( 26 ) 



with <7 x | y | 2 now skew Hermitian and the same scalar map as before. The superscripts W,( 2 ) and ^ denote 
operators acting on the Hilbert spaces of the first qubit, second qubit and the ancillary qubit respectively. 
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Fig. 3. (Left) Schematic of the ancillary quantum controller, that is allowed to interact with the 2-qubit system and the thermal 
bath via a tunable interaction. (Right) Prospective candidate for experimental implementation of the ancillary quantum controller 
t24l . Courtesy of Y. Nakamura, RIKEN, NEC, Japan. 



It is now seen that, 

[K X \ 2 , Kj] =c.J2 Vyil(9kb{ + g* k b k )\Z) (from eq. ©) 

k 

belongs to the control algebra generated by the additional vector fields introduced by the ancillary system, 
(i.e), [Ki, Kj] belongs to the Lie algebra generated by the control vector fields K±, ■ ■ ■ , Kg of the above 
system. Hence restructuring the system such that the linear span of the control vector fields and Lie 
algebra of the control vector fields coincide would ensure that the necessary and sufficient conditions 
given by equations d23l ) and (1241) are satisfied. 

VII. The Restructured Quantum Control System 

The ancillary qubit is primarily used to generate vector fields that can help decouple the system from 
the vector field Kj. The additional vector fields that can be generated with the help of ancillary system, 
help satisfy the necessary and sufficient conditions for decouplability as shown in this section. Let, 

j=l k 

denote the Hamiltonians of, qubits 1&2, environment and the ancillary system, 

tf^ElE^I (9kbl + g* k b k ) 
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the System + Environment decoherence Hamiltonian, 

H X = a« , H 2 = of , H 3 = of , H 4 = of 
the control Hamiltonians acting on qubits 1 and 2, 

H 5 = 2f 6 = of,H 7 = J,ofo^\H, = J 2 ofof 
the control Hamiltonians for the ancillary system, the Ising type coupling to qubits 1 and 2, and 

H 9 = ^o^(w k bi + wib k ) 

k 

the controllable interaction of the ancillary qubit with the environment. 

These 9-control Hamiltonians(f/i • • • #9) along with Hg and Hse govern the evolution of the system. 
The controls are implemented by the actual hardware(ancillary quantum controller) and corresponding 
fields (u\,---ug) with uj and iig being the strength of the Ising coupling. With this setup it is now 
possible to generate additional control vector fields by suitably manipulating the field strengths. With the 
additional controls generated, it is possible to satisfy the necessary and sufficient conditions (l23l l24l ). This 
way we are able to come up with the "restructured" quantum control system via the following "control 
pulse" maneuvers. For example, consider the following maneuver, with uq and ug, 

uq{t) = 1, and ug(r) = 0, for r G [0,i] 
, u&{ T ) = 0, and ug(r) = 1, for r € [t,2t] 
ue{r) = —1, and ug(r) = 0, for r € [2t, 3t] 
uq(t) = 0, and ug(r) = —1, for t € [3i, 4t] 

The corresponding unitary time evolution operator at the end of time instant At is given by, 

TJ(A t \ - e (-iHet) e (-iH g t) e (iH 6 t) e (iH t) 

= exp(-i[H 6 ,Hg]t 2 + 0(t 3 )) 

the series expansion by Campbell-Baker-Hausdorff formula. In the limit that t = dt — > 0. The effective 
direction of evolution is given by the commutator of the corresponding Hamiltonians, but to the second 
order in time. Hence we can devise a control vector field in the direction given by the commutators of 
the corresponding Hamiltonians Hq and Hg, where, 

[H 6 , H 9 ] = c.4 fc ) Yt w * b k + w th) 
k 
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where c is a real constant for a skew Hermitian Hq and Hg. In fact it is possible to generate a direction of 
evolution with arbitrary strength corresponding to repeated commutators of the Hamiltonians Hi-- - H$ 
of the physical system ( |26l ). The commutators of tensor product of operators are calculated according to, 

[A <g> B, C <g> D] = CA ® [B, D] + [A, C] ® BD 

With the control field H& we can generate the following direction in conjunction with the previous 
maneuver [HgjH^] = d ' Ji&z' ^ Oy and also, 

[[H 8 ,H 5 ], [H 6 ,H 9 ]} = Cl .[J 2 a^af\<j^ ^2(w k b\ + w* k b k )] 

k 

= c.aW<j^^(w k bl + w* k b k ) (27) 

k 

A similar maneuver between controls u^,uq and u%, generates the following direction of evolution, 

[H,,H,] = [af),J 2 af)a^]=c.a^af) (28) 

where c is a real constant for a skew Hermitian H±,H%. Again, from operating on equations (|27T ) and 
(EH)) we get, 

[[H 4 ,H 8 ], [[H 8 ,H 5 ], [H 6 ,H 9 ]\] = Cl [a^4 b \a^a^ £>4 + <h)\ 

k 

= ci\af\af\{af)\Y J ^kb\ + W lh) 

k 

= c. ff ( 2 )#.^( Ji ;4 + ^) (29) 
fc 

where 1^ is the identity operator on the ancillary subsystem. Hence we have generated an effective 

(2) 

controllable coupling between a y and the environment with the help of the ancillary qubit. It is important 
to note that the Hamiltonian so obtained by the above control maneuver now acts trivially on the Hilbert 
space of the ancillary qubit, a property which is found to be extremely useful. It is also possible to 

(2) 

generate the a x counterpart of the above coupling by a similar maneuver, given by, 

c.a^l( h lJ2( w kbl + w* k b k ) (30) 

k 

Again by a symmetric and identical argument we can generate a coupling between the environment and 
qubit 1, which is given by, 

c.aWjW. + w* k b k ) and c^J^. J2^ b l + w V>k) ( 31 ) 

k k 
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It can be seen that the above vector fields are what are required in equation (l25l) for the Lie bracket 
[iifm, to be contained within A + G. Now, noting that the constants c in the above equations can 
be controlled independently and arbitrarily, we can write the preliminary form of the actual control 
system which achieves disturbance decoupling. Gathering terms d29l)-(l3TT). we construct the following 
control system for given by equation (l32l ). In the following control system, the environment is 



approximated to be of single mode and of three energy levels II231 . The vast majority of the interaction 
energy is stored in the fundamental mode and first few energy states of the oscillator. The 24 restructured 



9\m) 



2 \ 2 



. E T a z } + 5>* 6 l** i£(*)> + E^ 6+ +9*b)\m) 

j=l k J j=l 

2 2 

+5> 1 ;<4 1) M> t + w*bym)) + E^4 1} ( w6t + w * b r\m) 

i=0 i=0 
2 2 

+ ^u 3l aM( W tf + w%y\m) + E^4 2) M + w*by\m) 

i=0 i=0 
2 2 



i=0 i=0 
2 2 

+ YjWPaPiwtf + w%ym) + ^u Si 4 1) 4 2) (^ + w * b ) l \^)) (32) 



i=0 i=0 



controls uiq ■ ■ ■ ni2,ii20 ' ' ' ^22 , ■ ■ • , ^30 ; ■ ■ ' ^82 could be thought of as "software" generated controls by 
manipulation of strengths of actual fields from system d26l ). Though, in the system above, we consider 
only first 3 states of the environment, it is possible to consider any finite number of environmental 
interaction terms by simply including additional control terms in system (l32l . via the same analysis. 
We are now in a position to use the new controls to decouple the scalar map from the environmental 
interaction. This restructured system satisfies the necessary and sufficient condition for decouplability 
because the disturbance vector field, Kj = Hse\Q, is contained within ker(dy) where y = (£|01)(10|£). 
To see this, one can evaluate L^y(t) and notice that it vanishes. The sufficient condition can be seen 
from the fact that [Kj,Ki] € G, where G is span of control vector fields of system (l32i Hence the 
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necessary and sufficient conditions, 

(i) Ki £ Ac ker(dy) 

(ii) [K I: Ki]eA + G, where G = span{#i • • • , K 2 a} (33) 

are satisfied. 

Scalability: It can seen that in order to perform the restructuring of a quantum system with finite number 
of qubits, it is necessary to employ only one ancillary quantum controller. This guarantees the scalability 
of the construction and decouplability properties of resulting system. 

In summary, the finite system and environment approximation has enabled us to come up with the 
control system whose coherence can be perfectly decoupled from the environmental interaction as shown 
in section (IVIIII ). By ensuring the maximum rank of control matrix j3, the controllability properties of 
the original 2-qubit quantum system on the manifold is fully preserved. In summary, given the existence 
of an invariant subspace A C ker(dy), the above conditions for decouplability can be summarized as, 



Open loop Uncontrolled 


[A, K ] c A 


[A, Ki] c A 


Controlled 


[A, Ko) CA + G 


[A, Ki] CA + G 



The analyses performed thus far also allows us to summarize the results in the form of the following 
corollary. 

Corollary VII.l. For any finite N-qubit open quantum system acted upon by arbitrary user generated 
control via Pauli matrices and under the influence of decoherence interaction Hse, the coherence between 
the basis states \i) and \ j) cannot be rendered immune to Hse without the action of the quantum controller 
which is the ancillary quantum system. 

VIII. Results 

The above system was simulated with 2, two-level interacting systems of interest and 1 ancillary qubit 
interacting with the environment. The goal was to study the effect of decoherence Hamiltonian on the 
coherence between |01) and 1 10) under arbitrary control as stated in the problem statement. We present 
the simulation results based on the above control strategy for two different strengths of decoherence 
interaction (red, for interaction strength 10 and black for no interaction) and different(both constant and 
time varying) control strengths. We approximate the environment by a single mode and first three levels 
for simplicity, as the vast majority of the environmental energy, modeled as infinite harmonic oscillators 
in different modes, is stored within the lowest mode and in lower energy states. This can be seen by the 
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Fig. 4. Open loop behavior of the 2 qubit system interacting with 3 level environment 



nature of the decoherence interaction Q and the coherent state of the harmonic oscillator |[25l . However, 
it can also be seen that, in order to include the higher order environmental interactions, one merely 
needs to add more control terms, following the same analysis. In this work, we consider the particular 
model of environmental interaction only to demonstrate the applicability of the technique and clarity of 
presentation. The initial coherence of the state between 1 10) and |01) is set to 0.5. The absolute value 
of the coherence with and without the decoherence interaction is presented along with the norm of the 
state in both the cases. Figure (@]) is the open loop behavior of the 2-qubit system and Figure (f5]) is the 
behavior of the decouplable system (l32l ). acted upon by the ancillary quantum controller. For the open 
loop system the coherence between states 1 10) and |01) is influenced in the presence of decoherence 
interaction. Though the coherence is irretrievably lost in reality, the periodicity is due to the finiteness 
of the system and environment. With the action of ancillary quantum controller and synthesis of control 
parameters as outlined in section [X] we see identical behavior of the coherence for different values of 
decoherence strength and useful control, thus effectively rendering the coherence between states 1 10) and 
1 01) immune to Hse- The behavior is seen to match for any set of analytical control functions, thus 
achieving perfect decoupling. In addition, the value of H in the Schrodinger equation was set to T, 
which gives rise to scaled time in the simulation. The slight deviation of the norm from ' 1' is an artifact 
of the numerical imperfections. 
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Fig. 5. Coherence between the basis states |01) and 1 10) of a 2-qubit system under the action of the ancillary quantum 
controller is preserved under arbitrary controls. The controls were chosen to be either constant, time varying or both. 
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IX. Quantum Internal Model Principle 

In this section we present the need for the Quantum Internal Model principle. Classical internal model 
principle for linear systems |[26l E71 outlines conditions for robust tracking in the presence of disturbance 
d(t). Consider the linear system with disturbance, 

x = Ax + Bu + E d d, (34) 

y = Cx + Du + F d d. (35) 

with the tracking error e = Cx + Du + F d d — r. The exosystem consists of the reference input and the 
plant noise, both generated by linear autonomous differential equations, 

r = A lr r,r(0) = r ; (36) 

d = A ld d, d(0) = d (37) 

with arbitrary initial states. The robust output regulation, where the tracking error is driven to zero, not 
only requires a dynamic state feedback but also that the controller mimic the exosystem in terms of its 
characteristic polynomial. Classical disturbance decoupling lfl"8l |fl9l GUI , on the other hand requires only 
the knowledge of system parameters and not the model of disturbance in order to completely decouple the 
output. However, quantum decoherence control, which is similar in formulation to classical disturbance 
decoupling is possible only with the knowledge of the environmental interaction, which is analogous to 
classical robust output regulation. This allows us to propose the quantum internal model principle with 
the following characteristics, 
Quantum and Classical Internal Model Principle 

• Quantum Internal Model principle aims at disturbance rejection with the help of a ancillary quantum 
controller and the knowledge of the model of interaction with the environment. 

• Classical Internal Model principle aims at perfect trajectory tracking via feedback, which involves 
the knowledge of the disturbance generator (as well as the desired trajectory generator) viz. the 
exosystem, 

Quantum and Classical Disturbance Rejection 

• Quantum disturbance decoupling, which is the underlying motivation of Quantum Internal Model 
requires complete knowledge of the model of the environment as well the corresponding model of 
decoherence, in the combined system+environment state space. 

• Classical disturbance decoupling only requires the model of interaction within the system's state 
space. 
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Hence the knowledge of environmental interaction with the system for complete decoupling makes 
quantum internal model principle salient and important within the framework of systems and control. The 
original 2-qubit system had to be augmented with the (internal) model of the environment entering via 
the control ug in equation d26l ) so as to restructure the vector fields to act non-trivially on the environment 
Hilbert space. Hence the knowledge of the model of interaction with the environment, i.e, the decoherence 
Hamiltonian Hse is essential to successfully controlling decoherence. Figure © depicts the nature of 



Restructured system with the 


(Internal) model of the 


interaction Hsb 


Original Open 




QuantumSystem 



Fig. 6. The original open quantum system acts as the skeletal structure for the larger restructured system. 

the original and restructured systems in the sense that the latter is larger and is derived from the original 
by taking into account the model of the environmental interaction Hse, whose closest classical analog 
is the disturbance generator denoted by A14. The Figure $7} outlines the schematic of control system for 
the decoupling problem, where the coherence measure for the controlled open quantum system and the 
corresponding closed system are identical. In summary, the structure of the system needed to be altered 
in order to, 

• Artificially induce coupling between qubits 1, 2 and the environment with the help of the ancillary 
qubit. 

• Generate vector fields in the higher order of the environment operator via a fast-action open-loop 
control. 

Hence it was necessary to modify the core system in more ways than one in order to perform decoupling. 
It is to be noted the above control strategy is a hybrid of fast-action open loop control and smooth analytic 
control (a and /?) in order to achieve perfect decoherence elimination. 

Conclusion 

In this article we visited an Internal Model Principle that is uniquely related to quantum systems in light 
of disturbance decoupling and decoherence control. The tensorial model of interaction of the quantum 
system with the environment can be skillfully exploited to completely decouple the system from the same. 
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Fig. 7. The difference between coherence measures from the open quantum system and the closed quantum systems must 
vanish. 



Such a result and its implication are first of its kind in the literature to the best of authors' knowledge. The 
ideas here presented could not only help further decoherence control but also influence the design of future 
quantum and classical control systems. In addition a framework for enhanced disturbance decoupling was 
laid wherein the entire control algebra can be used to effectively decouple a larger class of systems than 
just the linear span of the control vector fields. 



X. Appendix 

By following the proof of Theorem dIV.31) above as outlined in |[23l , we synthesize the control 
parameters a(£) and /3(£) for the System (l32l ). Let the restructured control vector fields in the system 
Hilbert Space be given by, 
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along with three environmental operators (I,D, D 2 ), we get 24 control vector fields, 



rW ^(1) a® a W 



a 



(2) 



{gi, ■■■92a} = M ',v y 

where the multiplication is carried out in the usual order, where D = (wb^ + w*b), is the displacement 
operator for Quantum Harmonic Oscillator. Define G = {<?i, • • • ,324}- The system and environmental 
identity operators are suppressed for ease of notations and they are assumed to be present where it is 
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clear from the context. The above control vector fields are generated by "software" in that, the control 
vector fields are produced by maneuvering action of fast action control pulses applied to the ancillary 
qubit enhanced open loop system (|26l ). 

In order to construct the invariant subspace for the restructured quantum control system with the control 
vector fields as above, it can be seen that the vector fields, 



Si 


= (4 J) 


+ 4 2) )I0, 




So 


= 4 1} c 


^4 2) ie> 




Ss 


= iie>, 






Si 


= (4 1} 


®4 2) -4^ 


5 4 2) )io 


s 5 


= (4 X) 




5 4 2) )I0 



commute with the vector field generated by the coherence operator, C|£} = (4 (&ax +4 ®<r^ ) |£) = 
|01)(10|£), which implies that L$ z y(t) = 0, or that the vector fields 6i are within ker(dy). It can also be 
seen that the the corresponding Lie brackets [Si, Sj] lie within A = span{<5i, • • • , S§}. 

Hence the invariant subspace for the above quantum system is identified to be generated by 5 Hermitian 
operators in the system's Hilbert space(but not all linearly independent for all the values of the states). 
This along with three commuting environmental operators(I, D, D 2 ) produces 15 vectors on the analytical 
manifold which span the invariant subspace for the restructured system. It can also be seen that since the 
5 system Hamiltonians do not always generate linearly independent vectors(5i, • • • ,65), the rank of the 
invariant subspace is dependent on the point £ and hence is singular. The methodology outlined in the 
proof to construct the control parameters a(£) and locally around the point £ works for non-singular 
invariant distribution A and non-singular control distribution G as well. We can now complete the basis 
for the tangent space T^(M) with the three commuting vector fields to A which do not belong to ker(dy), 



dl 


= (a« 


-4 2) )l£> 






= (a« 


(8 4 2) + 4 1} ° 


9^)10 


d. 


= (4 X) 




9 4 2) )l£> 



The commutation relations are as follows, 

•[SuS^eA; •[5 i ,g j ]eA + G 
• [Si,dj] G A, dj £ ker(dy); •[d i ,gj] G G 

Setup: 
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• Let K = rank{5i, • • • , J5} and {61, • • • , 6k} the corresponding vector fields with A = span{<5i, • • • , 5k}- 

• Let q be the minimum number such that rank{A, d\, ■ ■ ■ , d%} = rank{A, d\, ■ ■ ■ , d q }, q € 1, 2, 3 and 
let {di, • • • , d q } be the corresponding linearly independent vector fields with V q = spanjdx, • • • , d q }. 

• Let r be the minimum number such that rank{A, V q ,g±, • • • , 524} = rank{A, V q , g\, ■ ■ ■ , g r } with 
V = span{A,Vg,c/i,-- - ,g r } 

Let the vectors v 1 , • • • , v r be the linearly independent vectors of V according to the construction above. 
The Algorithm: 

• Solve the equation, 

24 K r 

^ 9jPji = ^ C ik V k + V K +i + ^ Ci k VK+q+k 
j=l k=l k=K+q+l 

for i = 1, ■ ■ ■ ,g with real coefficients This is obtained by rewriting the above equation as, 

[gi,--- ,924,-Vl, V K ,-V K +q+l, - ■ ■ ,~V r ]x 

[PlU " " " ; &"2AiiC%\i •■• ■, CiK, CK+q+1, • " " , Cir] T = ^X+i 

A least square solution to the above equation yields the local numerical values of the parameter 
f3ji £ R for all rows and first i = 1, • • • g columns. 

• Next, solve the equation Ylj=i 9jPji = Sfe=i c ifc^fc> for i = q + 1, • • • , 24 to obtain the rest of the 
parameters j3j{. This is again obtained by setting, 

[ffl, • • • ,524, -Vl, V r ] X {fin, ■ ■ ■ ,/3 2 4i,Cil, ' ' ' , C ir ] T = 

Hence the null space of the matrix [G, V] provides values for the parameters. 

• Finally the control parameters a are obtained by the solution to the equation, 

24 K r 

y~] ujgj + K = ^ c kVk + ^2 CkVk 

j=l k=l k=K+q+l 

or the least square solution to the matrix vector equation, 

[gi,-' ,924,-Vl, V K ,-VK+q+ir ■ ■ , ~V r ] X 

[di,-" , "24, Cl, • • • , CK, CK+l+l, • • • , Cr] T = —Kq 
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